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Abstract 

The  purpose  of  this  paper  is  to  minimize  the  energy 
consumption  for  a  single  link  communication  in  wireless 
sensor  network.  The  quadrature  amplitude  modulation 
(QAM)  and  auto  repeat  request  (ARQ)  with  finite  retry 
limit  and  various  packet  size  are  considered  in  the  system 
model.  We  formulate  the  optimization  problem  to  be 
a  nonlinear  integer  programming  with  modulation  order 
and  packet  size  as  variables.  To  solve  this  problem,  the 
continuous  relaxation  based  branch  and  bound  method  and 
the  successive  quadratic  programming  (SQP)  algorithm  is 
used.  Further,  a  variable  transformation  method  is  adopted 
to  convexificate  the  problem  and  guarantee  the  global 
optimization.  A  numerical  example  is  given  to  illustrate 
the  solution,  which  suggests  that  the  joint  optimization  of 
packet  size  and  modulation  order  for  various  communica¬ 
tion  distance  can  make  the  sensor  nodes  work  in  a  more 
energy  efficient  way. 


I.  Introduction 

In  recent  years,  wireless  networks  have  emerged  as 
an  efficient  real  time  communication  standard  for  both 
commercial  and  military  networks.  There  has  been  an 
urgent  need,  due  to  the  demand  for  applications  that  utilize 
data  from  multiple  remote  sources  to  classify  targets 
and  communicate  with  high  priority  responders,  to  make 
communications  on  wireless  networks  as  seamless  and 
uncorrupted  as  possible.  Unfortunately,  the  task  of  opti¬ 
mizing  communications  becomes  difficult  when  the  senor 
nodes  are  small,  and  have  a  low-battery  life  [Akyildiz  etal., 
2002],  Mission  critical  communications  in  the  battlefield 
could  be  disrupted  if  a  sensor  network  stops  working 
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suddenly.  In  situations  where  the  soldier  might  have  to 
identify  friend  from  foes,  intruders,  or  even  simply  obtain 
readings  from  sensor  nodes  as  to  the  conditions  of  a  par¬ 
ticular  unmanned  environment,  sensors  on  a  network  such 
as  distributed  Unattended  Ground  Sensors  that  have  failed 
or  are  failing  tend  to  contribute  erroneous  data  thereby 
providing  compromised  situational  awareness.  Thus,  a 
false  alarm  for  a  target-rich  zone  could  be  returned  to 
the  warfighter,  who  would  enter  the  zone,  unprepared  for 
conflict.  With  these  concerns  in  mind,  research  into  the 
realm  of  energy-efficient  communication  strategies  from 
the  perspective  of  cross  communication  layer  optimization 
has  been  conducted.  Many  energy-efficient  strategies  have 
been  proposed  in  different  layers  or  by  the  use  of  cross¬ 
layer  design,  which  is  well  surveyed  in  [Goldsmith  etal., 
2002], 

In  physical  layer,  the  optimization  of  modulation 
order  is  proposed  in  [Cui  etal.,  2005]  to  achieve  a  balance 
between  circuit  energy  consumption  and  transmission  en¬ 
ergy  consumption.  The  work  is  extended  to  the  multiple 
access  scenario  by  [Cui  etal.,  2004].  Further,  [Prabhu  etal., 
2007]  minimized  energy  consumption  by  joint  modulation 
and  power  adaptation  in  consideration  of  auto  retransmis¬ 
sion  request  (ARQ)  in  data  link  layer.  It  assumed  a  simple 
ARQ  scheme  with  infinite  retransmission  limit  to  simplify 
the  problem.  Unfortunately,  considering  only  optimization 
of  the  modulation  order  omits  features  of  the  data  itself 
such  as  cyclic  redundancy  codes  and  perhaps  more  com¬ 
plex  strategies  such  as  data  compression,  which  might 
allow  for  high  fidelity  transmission  at  a  lower  energy 
cost.In  data  link  layer,  energy  efficiency  based  packet- 
size  optimization  is  implemented  in  [Sankarasubramaniam 
etal.,  2003],  Once  again  packet  size  optimization  tech¬ 
niques  based  on  parameters  such  as  channel  conditions 
does  not  provide  a  very  complete  energy  efficient  solution 
since  without  considering  the  setup  of  the  under-lying 
hardware,  there  is  a  fundamental  limit  to  how  energy- 
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including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

01  DEC  2008 

2.  REPORT  TYPE 

N/A 

3.  DATES  COVERED 

4.  TITLE  AND  SUBTITLE 

Optimization  of  Energy  Consumptions  in  Sensor  Networks 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Department  of  ECE  Stevens  Institute  of  Technology  Hoboken,  NJ  07030 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

See  also  ADM002187.  Proceedings  of  the  Army  Science  Conference  (26th)  Held  in  Orlando,  Florida  on  1-4 
December  2008,  The  original  document  contains  color  images. 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

uu 

18.  NUMBER 
OF  PAGES 

6 

19a.  NAME  OF 
RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


And  the  average  transmission  times  for  each  packet 


efficient  transmissions  can  get,  constrained  by  the  physical 
configuration  of  the  system.  By  introducing  the  cross¬ 
layer  design  scheme,  [Kwon  etal.,  2006]  considered  power 
control  in  physical  layer,  retransmission  limit  in  datalink 
layer  and  routing  protocol  in  network  layer  to  maximize 
the  lifetime  of  WSN.  However,  the  modulation  order  and 
packet  size  are  not  involved  in  the  optimization.  A  general 
optimization  model  for  bottom  three  layers  of  WSN  is 
outlined  in  [Cui  etal.,  2007], 

In  this  paper,  we  focus  on  the  cross  layer  energy 
minimization  for  a  single  link  of  WSN  in  consideration  of 
the  modulation  order  in  physical  layer  and  packet  size  in 
data  link  layer.  An  ARQ  scheme  with  finite  retransmission 
limit  is  assumed.  And  the  transmission  power  is  adjusted 
according  to  the  modulation  scheme  and  packet  size  to 
approach  the  target  packet  loss  rate  with  the  given  retrans¬ 
mission  limit.  This  work  is  a  basis  of  the  whole  network 
energy  minimization  and  has  so  far  been  overlooked. 
The  contribution  of  this  paper  is  to  establish  the  energy 
consumption  model  for  a  single  WSN  link  with  MQAM 
modulation  and  finite  retransmission  ARQ  scheme  and 
jointly  optimize  the  modulation  order  and  packet  size.  As 
a  typical  nonlinear  integer  programming,  the  optimization 
problem  can  be  solved  by  the  continuous-relaxation-based 
branch-and-bound  method  [Li  etal.,  2006].  However,  it  can 
be  found  that  the  problem  is  non-convex  even  after  con¬ 
tinuous  relaxation.  Variable  transformation  method  is  used 
to  convexificate  the  continuous  problem  [Li  etal.,  2005], 
Then  the  global  optimization  can  be  obtained  by  general 
nonlinear  convex  optimization  algorithms,  such  as  succes¬ 
sive  quadratic  programming  (SQP)  method  [Bazaraa  etal., 
2006], 

The  rest  of  the  paper  is  organized  as  follows.  Section 
II  describes  the  system  model  and  section  III  formulats  the 
optimization  problem  as  a  nonlinear  integer  programming 
and  discusses  the  optimization  algorithm  and  convexifica- 
tion  method.  Some  numerical  results  are  given  in  Section 
IV  and  the  conclusion  is  drawn  in  section  V. 


II.  System  model 

A  single  communication  link  connecting  two  wireless 
sensor  nodes  is  considered.  Supposed  that  the  information 
bits  are  divided  into  packets  with  L  =  Lh+Lp  bits/packet, 
where  Lh,  and  Lp  are  the  fixed  overhead  bit  number  and 
the  payload,  respectively.  With  a  simple  ARQ  protocol, 
each  packet  will  be  transmitted  repetitively  until  it  is  re¬ 
ceived  successfully  or  the  maximum  transmission  number 
Nr  occurs.  If  the  target  packet  loss  rate  is  pioss,  the  packet 
error  rate  (PER)  of  each  transmission  can  be  calculated  by 
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For  each  packet  transmission  in  the  link,  the  nodes 
will  spend  time  Ton  in  the  active  mode  for  transmission  or 
reception,  then  return  sleep  mode  in  a  transient  duration  of 
Ttr.  Assumed  that  the  energy  consumption  in  sleep  mode 
is  0,  the  total  energy  consumption  E  required  to  send  a 
packet  successfully  can  be  calculated  by 

E  =  N{PonTon  +  PtrTtr) 

=  N[(Pt  +  Pc0)Ton  +  PtrTtr],  (3) 

where  Pon,  Ptr  are  power  consumption  values  for  the 
active  mode  and  the  transient  mode,  respectively.  Pt  and 
Pc o  are  the  transmission  power  and  circuit  power  in  active 
mode,  respectively.  Specially,  Pc0,  which  consists  of  two 
parts  of  power  consumptions,  can  be  given  by 

PcO  Pamp  T"  Pc 

=  aPt  +  Pc,  (4) 

where  Pamp  and  Pc  are  the  power  consumption  of  radio 
frequency  (RF)  power  amplifier  (PA)  and  other  circuits, 
respectively,  a  =  |  —  1  with  r]  the  drain  efficiency  of  the 
RF  PA  and  £  the  signal  peak-to-average  ratio  (PAR). 

Based  on  the  model  outlined  above,  the  energy 
consumption  for  the  successful  communication  of  each 
information  bit  in  the  link  can  be  given  by 

E  N[(^jPt  +  Pc)Ton  +  PtrTtr\ 

PpiX  Pioss)  (P  Ph){\  Pioss) 

(5) 

In  physical  layer,  we  assume  that  the  bandwidth  is  B 
and  the  data  rate  is  77,  bits/symbol.  It  is  obvious  that  the 
packet  size  L  should  be  an  integer  times  of  77,  to  get  more 
energy  efficient.  Thus  the  symbol  number  in  each  packet 
is  Ls  =  L/ti,  and  the  time  duration  of  a  packet  is 

Ton  =  LSTS  «  (6) 

JD 

where  Ts  is  the  time  duration  of  a  symbol.  And  the  target 
symbol  error  rate  (SER)  can  be  calculated  by 

Pser  =  1  -  (1  -  pper )t* 

=  1-(1  -pfas)^-  (?) 

To  simplify  the  model,  forward  error  control  (FEC) 
is  not  considered  in  our  work.  The  information  bits  are 


Pper  =  P^ss- 


MQAM  modulated  with  rectangle  constellation  of  size 
M  =  2rb.  The  PAR  of  the  modulated  symbol  is 


2rb>2  -  1 

C  =  Q _ 

^  2rb/2  +  1  ’ 


(8) 


Supposed  that  the  system  works  on  an  additive  white 
Gaussian  noise  (AWGN)  channel.  The  bound  of  the  SER 
of  MQAM  is  [Proakis,  2001] 

PSER  <  4 Q{\j x) 

<  (9) 

where  Q(x)  =  fx(  1/ \/2n)e~u  l2du\  7  is  the  signal  to 
noise  ratio  (SNR)  with  7  =  Pr/(2Ba2Nf),  where  Pr  is 
the  received  signal  power,  cr2  is  the  power  spectral  density 
of  the  AWGN,  and  Nf  is  the  receiver  noise  figure  defined 
as  Nf  =  Ntotai/ (2f?<r2)  in  which  Ntotai  is  the  power  of 
the  noise  introduced  by  the  receiver  front  end. 

Considering  a  propagation  path-loss  model,  the  trans¬ 
mission  power  can  be  calculated  by 


Pt  =  GdPr ,  (10) 


III.  Problem  Formulation  and  Resolution 
A.  Problem  Formulation 

The  objective  of  this  paper  is  to  minimize  the  energy 
consumption  for  the  communication  of  each  bit  in  a 
single  WSN  link  considering  the  constraint  of  maximum 
transmission  power  and  average  delay.  The  problem  can 
be  written  as 


min  Eb,  (16) 

Tb,Ls 

S.t.  Pt  <  Pmaxi 

E  —  Pnax  5 

LsVb 

n >  ls  e  z+, 

where  Pmax  and  rmax  are  the  maximum  transmission 
power  and  maximum  average  delay,  respectively;  r  is  the 
average  delay  for  a  packet  to  be  received  successfully, 
which  can  be  approximately  calculated  by 

T  =  N(^+Ttr).  (17) 

JD 


where  Gd  =  G\dnMi  is  the  power  gain  factor  in  which 
Mi  is  the  link  margin  compesating  the  hardware  process 
variations  and  other  additive  background  noise  or  inter¬ 
ference  and  Gi  is  the  gain  factor  at  d  =  lm,  which  is 
dependent  on  the  antenna  gain,  carrier  frequency  and  other 
system  parameters;  n  is  the  path-loss  factor. 

According  to  (9)  and  (10),  by  approximating  the 
bound  as  an  equality,  the  required  transmission  power  is 

Pt~  ^GdNf Ba2(2rb  .  (11) 

•J  PSER 

By  taking  (6),  (8)  and  (11)  into  (5),  we  can  get 


Therefore,  the  delay  constraint  corresponds  to  the 
upper  bound  on  Ls  with  Ls  <  Lsmax  =  [B(rmaJ|~TtT')], 
where  [_x\  denotes  the  maximum  integer  less  than  or  equal 
to  x. 

Define  x  =  [27,  X2]  with  x\  =  Tb  and  27  =  Ls,  the 
optimization  problem  can  be  easily  formulated  to 

min  Eb  =  /0(x),  (18) 

X 

s-t.  /i(x)  <  0, 

f'2  (n)  <  0, 

lb  <  x  <  ub, 

x  e  z2+, 


Eb  = 


Ls 


Ls^b  L} 


-{Ai(2Tf  —  l)2  In  - 


1-(1  -pLsY 


x4_3  , 

+^2+/}, 


(12) 


with 


AGdNf<r2 

(13) 

VlX-PlJaa) 

Pc 

A2—  1  , 

(14) 

B{1-pLs) 

PrrTtr 

A3=  1  • 

(15) 

1  -Plols 

with 


/o(x)= 


Ai(2  21  -  l)2x2(i(x2)  +  A2x2  +  A3 
x\x2  -  Lh 


f1(K)=-GdNfBa2(2^  -  1  )h(x2)  -  Pn 


f'2  (x)  = 

X\X2 

h(x2)= In  — 


1 

Lh' 
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(19) 

(20) 
(21) 
(22) 


where  lb  =  [f'hrnin  •  and  ub  =  [t'hrnax  •  Lsrnax\  are 

the  lower  and  upper  bounds  of  x,  respectively.  Assuming 
that  4-QAM  is  the  smallest  constellation  size  available 
[Cui  etal.,  2005],  we  set  r6min  =  2,  rmax  =  +00, 


Lsmin  —  1  and  Lsrnax  [ 


B(Tmaa.-Ttr) 
N 


J- 
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B.  Continuous-based  Branch-and-Bound  method 


Obviously,  problem  (18)  is  a  nonlinear  integer  pro¬ 
gramming  (NIP),  which  generally  can  be  solved  by  the 
continuous-relaxation-based  branch-and-bound  method, 
whose  detailed  process  is  outlined  below  [Li  etal.,  2006], 

The  algorithm  starts  by  finding  an  optimum  solution 
x*  and  the  minimized  object  value  fu,  of  the  continuous 
relaxation  of  the  NIP  problem.  If  x*  is  integral,  then 
it  is  also  optimum  for  the  NIP  problem.  Otherwise,  we 
approximate  the  integer  optimum  solution  as  x*  =  x* 
with  corresponding  object  value  fub  =  /0(x*).  If  the 
relative  error  c  =  ubfJlb  caused  by  the  approximation  is 
less  than  a  given  tolerance  e,  x*  is  regarded  as  solution  of 
the  NIP  problem.  Otherwise,  select  a  fractional  variable 
of  x*,  e.g.  x*.  Two  new  subproblems  are  generated  by 
adding  variable  constraints  Xi  <  x*  and  Xi  >  x*  +  1, 
respectively.  For  each  subproblem,  we  apply  the  same 
method  to  find  fu,  and  fub.  Choosing  the  subproblem  that 
gives  the  minimum  fib,  we  can  continue  to  divide  it  until 
f  <  e  or  the  problem  is  infeasible. 

C.  Variable  transformation  convexification 

It  can  be  found  that  the  above  algorithm  is  based  on 
the  solution  of  the  continuous  relaxation  of  the  NIP  prob¬ 
lem.  However,  even  after  continuous  relaxation,  problem 
(18)  is  non-convex.  Therefore,  the  global  optimization  can 
not  be  guaranteed  by  convex  optimization  algorithms.  For¬ 
tunately,  we  find  that  this  problem  can  be  convexificated 
in  some  region  by  a  variable  transformation  defined  by 


y  =  t(x), 

(23) 

l 

yi  =  ti(xi)=xl , 

(24) 

_ ]L 

y2  =  t2(x2)=x2  4. 

(25) 

Obviously,  the  transformation  is  an  one-to-one  map¬ 
ping,  so  the  continuous-relaxation  of  problem  (18)  is 
equivalent  to  [Li  etal.,  2005] 

min  Eb  =  f0(  y),  (26) 

y 

s.t.  /l(y)  <  0, 

h(y)  <  o, 

lby  <  y  <  uby, 

with 

/o(y)  =  /o(t-1(y))  +  p-, 

L h 

=  Y^)[Al9l{yi)Hm)  +  $  +  (27) 


A(y)  =  /i(t  Ay)), 

=  ^GdN fBa2g2(yx)h{y2)  -  Pmax, 

(28) 

My)  =  /2(t_1(y))  =  z( y)  -  1, 

(29) 

,  ,  (2^-l)2 

9i(yi)  —  s  > 

(30) 

Vi 

32(2/1)  =  2yi  -  1, 

(31) 

~h{y2)  =  h(t2l(y2))  =  In 

(32) 

1"(1  -Piol.)v> 

z{y)  =  LfM, 

(33) 

i/i 


where  lby  and  uby  are  the  lower  and  upper  bounds  of 
y,  respectively. 

To  analyze  the  convexity  of  problem  (26),  two  impor¬ 
tant  properties  of  convex  functions  will  be  used:  (1)  The 
sum  of  convex  functions  are  convex;  (2)  The  product  of 
log-convex  functions  are  convex  [Bazaraa  etal.,  2006]. 

By  calculating  the  second-order  derivations  or  Hes¬ 
sian  matrices,  it  is  easy  to  prove  that  the  functions 
lngi(a;).  In  g2(x),  In  jM.,  In  and  z( y)  are  convex. 
In  addition.  In  y3_  is  nondecreasing,  therefore.  In 

is  convex.  Further,  the  calculation  of  :l  j£ —  indicates 
that  h(x)  is  convex  when  1  —  (1  —  pl/)SS')x  <  e-4, 

which  is  equivalent  to  Ls  >  L’smin  = 
in  our  model,  where  \x~\  denotes  the  minimum  integer 
larger  than  or  equal  to  x.  Thus,  according  to  the  property 
(1)  and  (2),  /0(y),  /i (y)  and  f2( y)  are  convex,  i.e. 
problem  (26)  is  convex,  in  the  region  of  [lby,uby],  with 
lby  =  [21/3,L~^/ai]  and  uby  =  [+oo, 

In  other  words,  the  continuous  relaxation  of  en¬ 
ergy  minimization  problem  (16)  can  be  transformed 
to  an  equivalent  convex  problem  (26)  when  Ls  £ 
[L'smivx  LSinax}-  So  conventional  nonlinear  convex  opti¬ 
mization  algorithms,  such  as  successive  quadratic  pro¬ 
gramming  (SQP)  method  [Bazaraa  etal.,  2006],  can  be 
adopted  to  get  the  continuous  optimum  point  and  the 
branch-and-bound  method  can  be  used  to  search  the 
integer  optimum  value  of  rb  and  Ls  in  this  region. 
Moreover,  L'smin  is  very  small  in  general,  e.g,  when 
Pioss  =  10-4,  L'smin  =  1  if  Nr  =  1, 2  and  L'smin  =  3 
if  Nr  =  3.  Therefore,  the  exhaustive  search  method 
can  be  used  to  minimize  the  energy  consumption  for 
Ls  £  [Lsmim  L'smin}.  Then  a  global  optimization  can  be 
obtained  by  comparing  the  minimization  results  in  the  two 
regions. 
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(b)  Transfomed:  E 5  v.s.  y2 
Fig.  1.  Objective  function  with  7*5  =  2  at  d  =  50m 

IV.  NUMERICAL  RESULTS 

For  a  numerical  example,  the  following  values  are 
configured  according  to  [Cui  etal.,  2005],  The  system 
bandwidth  B  =  10kHz,  circuit  power  Pc  =  210mw, 
transient  power  Ptr  =  lOOmw,  transient  time  Ttr  =  5/rs, 
G\  =  30dB,  path  loss  factor  n  =  3.5,  tr2  =  - 1 71dbm/Hz, 
Nf  =  lOdB,  r/  =  0.35.  The  overhead  bit  number  is  fixed 
to  Lh  =  32bits.  The  maximum  transmit  power  Pmax  = 
250mw  and  maximum  average  delay  rmaa;  =  500ms.  The 
target  packet  loss  rate  is  pioss  =  10-4  and  the  maximum 
retransmission  number  is  Nr  =  2. 

Fig.  1(a)  and  Fig.  1(b)  show  the  objective  functions  at 
d  =  50m  before  and  after  the  variable  transformation  with 
7-;,  =  2  and  variable  packet  size.  It  gives  us  an  example 
to  show  that  the  non-convex  objective  function  can  be 
convexificated  efficiently  by  the  variable  transformation. 
Fig.  2  and  Fig.  3  provides  the  optimum  values  for  the  mod¬ 
ulation  order  r b  and  packet  size  L  =  Lsrb,  respectively, 
with  different  communication  distance  d,  where  the  SQP 
algorithm  and  branch-and-bound  method  with  e  =  10-4 
is  used  to  solve  the  optimization  problem. 


Fig.  2.  Optimum  modulation  order  v.s.  Communication  distance 


Fig.  3.  Optimum  packet  size  v.s.  Communication  distance 


Fig.  4  compares  the  energy  consumption  of  our  joint 
optimization  system  with  that  of  the  system  with  fixed 
packet  size  and  optimized  modulation  order  as  proposed 
in  [Cui  etal.,  2005].  The  comparison  is  measured  by  the 
use  of  energy  saving  percent  defined  as 

P  =  Ef~  E°vt  x  100%,  (34) 

Ef 

where  Eopt  and  Ef  are  the  energy  consumptions  with 
joint  optimization  and  modulation  order  optimization  only, 
respectively. 

Observing  Fig.  4,  one  finds  that  the  joint  optimization 
saves  about  6%  energy  at  d  =  50 m  in  comparison  with 
the  fixed  packet  size  L  =  4000  and  save  about  4%  energy 
at  d  =  1  to  in  comparison  with  L  =  500.  The  experiment 
results  indicate  that  the  joint  optimization  can  make  the 
sensor  nodes  work  in  a  more  energy  efficient  way. 
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Fig.  4.  Energy  saving  by  joint  optimization  (in  comparison  with  fixed 
packet  size  and  optimized  modulation  order) 

V.  Conclusion 

In  this  paper,  we  derived  an  energy  consumption 
model  for  WSN  considering  the  parameters  in  the  physical 
and  data  link  layers.  Based  on  the  model,  a  NIP  problem 
is  formulated  with  the  modulation  order  and  packet  size  as 
variables  and  solved  by  the  use  of  continuous-relaxation 
based  branch-and-bound  method,  SQP  algorithm  and  vari¬ 
able  transformation  convexification.  It  was  found  that 
the  proposed  joint  optimization  can  improve  the  energy 
efficiency  of  WSN. 
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